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ABSTRACT 



Aims. We study the uniformity of the distribution of compact flat-spectrum AGN on the sky and the evolution of their relativistic jets with 
^ ■ cosmic epoch. 

tf^l Methods. A complete sample of compact extragalactic radio sources at 15 GHz was recently compiled to conduct the MOJAVE program. The 
C MOJAVE sample comprises 133 radio-loud flat-spectrum AGN with compact relativistic outflows detected at parsec scales. 

Results. Analysis of the population of flat-spectrum quasars of the sample reveals that the pc-scale jets of quasars have intrinsic luminosities 
in the range between ~ 10 24 WHz ' and ~ 10 27 WHz~' and Lorentz factors distributed between 3 S> y S> 30. We find that the apparent speed 
(or Lorentz factor) of jets evolves with redshift, increasing from z ~ to z ~ 1 and then falling at higher redshifts (z ~ 2.5) by a factor of 2.5. 
The evolution of apparent speeds does not affect significantly the evolution of the beamed luminosity function of quasars, which is most likely 
>sz , to be dependent on the evolution of radio luminosity. Furthermore, the beamed radio luminosity function suggests that the intrinsic luminosity 
J~i function of quasars has a double power-law form: it is flat at low luminosities and steep at high luminosities. There is a positive evolution of 
Oh quasars at low redshifts (z < 0.5) and strong negative evolution at redshifts > 1.7 with space density decline up to z ~ 2.5. This implies that the 
^ powerful jets were more populous at redshifts between 0.5 and 1.7. We show that the evolution of compact quasars is luminosity dependent 
5^ ' and it depends strongly on the speed of the jet suggesting that there are two distinct populations of quasars with slow and fast jets which evolve 
[ differently with redshift. 
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1 . Introduction and radio galaxies from larger angles. It was successfully used 

to describe a diverse appearance of observed char acteristics 

The advantage of studying the cosmological evolution of ra- f radio sources and their evolution Jli7r7& Padovanilfl99i 

dio sources free of dust obscuration was realized by|Longair| |j ac kson & WaiiFl999tlGrimes et al.l2004 at different radio fre- 

il966|). He demonstrated that the cosmic evolution of low que ncies. The orientation of the jet and its beamed emission are 

and high luminosity radio sources sampled at low radio fre- two important keys t0 identify that the BL Lacs and quasars 

quency (178 MHz) are different with later ones having much m bea med parent objects of FRI and FRII radio sources. At 

stronger evolution. The morphological difference between FRI low f reqU encies (178 MHz), these sources appear as extended 

and FRII radio sources (|Fanaroff & Riley||l974|) and their di- (i be-dominated) st eep-spectrum FRI and FRII radio sources, 

vision into low and high luminosity classes below and above Iwillott et alJ lEool used the combined sample of 3CRR, 7C 

ii78MHz ~ 10 25 WHz [ sr 1 was used by | WaU Pt alJ and 6CE surveys (see references in Willottetal. 2001) to model 

to show that the population of FRI radio sources does not the radio i um i nosity function of FRI/FRII radio sources. They 

evolve strongly with cosmic time whilst the evolution of FRII showed that a dual-population model of luminosity function 

sources is much stronger. The orientation-based unification fits the data well? but it requires a differential positive density 

scheme (|Scheuer||l987t |Barthel||l989|) suggests that a single evolution between z ~ and z ~ 2. No evidence of a negative 

parent population of radio sources appears as radio galaxies evolution (redshift 'cut-off') is found at high redshifts. 
and quasars when viewed at different angles to the line of 

sight, with quasars beeing seen from 45° to the line of sight The surveys at frequencies > 1.4 GHz include both lobe- 

and core-dominated sources in which the jet beaming effect 

Send offprint requests to: T. G. Arshakian appears to be dominant. These radio sources enhanced by 

* On leave from Byurakan Astrophysical Observatory, Byurakan Doppler beaming of the jet are identified with flat-spectrum 

378433, Armenia and Isaac Newton Institute of Chile, Armenian quasars and BL Lacs. At these frequencies, the cosmologi- 

Branch cal evolution of flat-spectrum sources is more complex, in- 
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volving the evolution of both non-beamed and beamed emis- 
sion. In addition to the luminosity/density evolution it also 
may depend on the evolution of relativistic effects with cos- 
mic time. The evol ution of flat-spectrum quasars was stud- 
ied in several works (IDu nlop & Pe acock! 1 990t [Jackson & Wall 
ll999tlCirasuolo et all2005t IWall et alJl2Q05l) . I.Tackson & Wall 
( 1999) showed that the evolution and beaming of powerful 
radio sources may be described by a dual-population uni- 
fication scheme for FRI and FRII sources which, depend- 
ing on the viewing angle, appear as lobe-dominated steep- 
spectrum and core-dominated flat-spectrum sources. Using 
data from a comb i ned sample of radio sources at 2.7 GHz 
iDunlop & Peacockl Jl99(]l) showed that both pure luminosity 
and luminosity/density evolution fit the observed redshift and 
source-count data. They presented clear evidence of a redshift 
cut-off for both steep- and flat-spectrum quasars at high red- 
shifts, z ~ 2. A high redshift cut-off is al so evident in the 
Parkes quarter-Jansky flat-spectrum quasars JWall et alJl200l 
and in the sample of 352 faint flat-/steep-spectrum quasars 
JCirasuolo et alJl2005l) - with radio fluxes >1 mJy at 1.4 GHz 
- selected from the radio (Faint Images of the Radio Sky at 
Twenty centimeters) and optical (2dF QSO Redshift Survey) 
surveys. They found an indication of a significant negative 
evolution for faint quasars at z > 1.8 consistent with nega- 
tive evolution of bright flat-spectrum radio quasars at 2.7 GHz, 
and inconsistent with no-evolution of steep-spectrum bright ra- 
dio quasars at 151 MHz JWillott et ai1ll998l) . The reduction of 
the comoving space density of quasars at high redshifts is also 
found in X-ray and optically selected samples (e.g., Schmidt 
et al. 1991; Fan et al. 2001 and references therein; Hasinger 
2004; Silverman et al. 2005) indicating that the redshift cut- 
off of optically-selected quasars is not a selection effect due to 
obscuration by dust. 

It is not clear whether the redshift cut-off of flat-spectrum 
sources is present in high radio frequency samples selected on 
the basis of beamed emission and whether the cosmological 
evolution of beamed sources depends on the properties of the 
jets such as the intrinsic luminosity and Lorentz factor of the jet 
(the analytical relationships be tween the latter parameters were 
studied bv lLister & Marscheil(ll997h for simulated flux-limited 
samples of relativistic jets). To examine these problems we use 
the (first) complete sample of 133 brig ht AGN at 15 GHz which 
has been compiled recently by iLister & Homan (2005) as an 
exte nsion of the Very Long Baseline Array (VL B A) 2 cm sur- 
vey llKellermann et al.ll998HZensus et alT20 02) and other pro- 
grams. This sample of flat-spectrum AGN was designed for a 
long-term observational program (named the MOJAVE survey) 
to study the structure and evolution of compact jets. The ra- 
dio emission of AGN originates in the compact jets at scales 
from a few parsecs to tens of parsecs. Although the sample 
is small, includes relatively bright sources and is restricted to 
Z ^ 2.7, it provides a wealth of information about radio emis- 
sion, structure and kinematics of parsec-scale jets. The com- 
pleteness of the sample makes it possible to investigate the 
cosmological evolution of flat-spectrum AGN in the context of 
evolution of the jet parameters. In Sect. 2, the complete sam- 
ple of flat-spectrum AGN is presented. In Sect. 3, we investi- 
gate the homogeneity of the MOJAVE sources on the sky. The 



source counts of AGN, evolution of apparent speeds of jets, 
and luminosity function and its evolution with cosmic epoch 
are examined in Sect. 4. Section 5 involves the summary of re- 
sults. In this paper we use a flat cosmology (D. m + £2a = 1) with 
non-zero lambda, £2 m = 0.3 and Ho = 70kms~' Mpc -1 , 

2. The MOJAVE sample 

The main selection criterion is the flux density limit at 15 GHz: 
all variable sources with measured VLBA flux densities ex- 
ceeding 1.5 Jy (2 Jy for southern sou rces) at any epoch since 

j ir — 

1994 are included in the sample (Lister & Homan 2005). It is 
restricted to declinations larger than -20° and to Galactic lat- 
itudes \b\ > 2.5°. There are 96 sources in common with the 
2 cm survey, which is not entirely complete down to the flux 
density limit. A sample of additional 37 AGN was assembled 
from other recent high-frequency surveys. The final sample 
comprises 133 radio sources, all active galactic nuclei: radio- 
loud, core-dominated and flat-spectrum 1 (see Lister & Homan 
2005, for more detailed description of the sample and its com- 
pleteness). The MOJAVE sample is based on the selection of 
highly beamed radio sources. The majority of the AGN have 
a core-jet structure on parsec-scales, and many of them show 
superluminal motion in their jets suggesting that the bulk of 
the beamed emission at 15 GHz originates in the relativistic jet. 
Some nearby radio galaxies with two sided jets are included 
in the sample because of bright (unboosted) radio emission. At 
present, 116 (87%) of the 133 AGN have spectroscopic red- 
shifts: 94 of 94 (100 %) quasars, 15 of 22 (68 %) BL Lacs and 
7 of 7 (100 %) radio galaxies. The remaining nine sources have 
no optical identification. The radio luminosities of compact 
AGN at 15 GHz (Z45) vary over (10 23 to 10 29 ) WHz 1 over 
the redshift range from to 2.7 (Fig. [3- The apparent speed 
of th e fastest component of t he jet rated as 'good' or 'excel- 
lent' jKellermann et al.ll2004l) is measured for 66 quasars, 15 
BL Lacs and six radio galaxies. The maximum total flux den- 
sities observed for each source during the VLBA monitoring 
campaign are used throughout the paper. 

3. Sky distribution 

If the sample is free from selection effects, then we should ex- 
pect that the radio sources are distributed uniformly in the sky. 
The distribution of 133 AGN on the celestial sphere is shown in 
Fig. [2] The distribution of radio sources looks patchy with few 
clusters and voids in the sky. To test the uniform distribution 
of AGN we m easure a two-p oint angular correlation function 
dPeeblesll98(iFnf a ntell994 



u>(8) = 



2N gg N T 



(1) 



where N g is the number of AGN in the MOJAVE survey, N gg is 
the number of distinct pairs of AGN with angular separations 
from 6 to + Ad, N T is the number of random points generated 
in the area of the survey, and N gr is the total number of random 
points in the same annuli 9 to + AO around AGN. Note, that 
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Fig. 1. The apparent luminosity at 15 GHz versus redshift for 
the 116 AGN with known redshifts. The filled circles, dia- 
monds and triangles are quasars, BL Lacs and galaxies, respec- 
tively. The solid and dotted lines show the flux density limits 
for northern (1.5 Jy) and southern (2 Jy) radio sources respec- 
tively, and the dotted-dashed line corresponds to the maximum 
flux density (41 Jy at a redshift 0. 158) observed in the MOJAVE 
sample. 
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Right Ascension 

Fig. 2. The distribution of 133 AGN in the region of the sky 
restricted by 6 > -20° (the dashed line). The open circles are 
sources with no optical counterpart, other symbols are the same 
as in Fig.[0 The dotted streams indicate the galactic plane ex- 
clusion strip, \b\ > 2.5°. 

only 4% of the sources {2nd I An, where d = 5°) are missed 
from the sample as a result of the Galactic plane exclusion. 
This has little effect on number statistics of AGN therefore we 
do not take it into account. The angular correlation function for 
133 AGN is shown in Fig. [3] The correlation function remains 
near a>(0) = within the limits of errors, which is indicative 
that the distribution of MOJAVE sources is uniform over the 
large angular scale from 5° to 40°. 

4. Space distribution: evolution of compact AGN 

At high radio frequencies the observed flux density is produced 
predominantly in the relativistic jet, 

L 



4nR 2 



(l+z> 



(2) 



where L is the beamed (apparent) luminosity, R is the comoving 
distance at the redshift z, and a is the spectral index (a = is 



Ans>le 9 



Fig. 3. The correlation function against angular separation for 
133 AGN. The uncertainty in u>(0) in each bin is given by the 
Poisson error, yfTTai(Uj/ JW^ (Infante 1994). 



assumed throughout the paper). The relativistic Doppler effect 
causes the rest-frame (intrinsic) luminosity (£) of the source 
to appear enhanced towards the observer, L = 5 P where p 
2 i > for a steady-state jet and 6 is the Doppler factor which is 
a function of the jet viewing angle (6) and the speed of the jet 
P = (1 - 1 1 y 2 ) 1 12 (in units of the speed of the light) and y is the 
Lorentz factor, 



6 = 



1 



y(l -p cos 0) 



(3) 



Therefore, the source counts and cosmological evolution of 
AGN from the MOJAVE sample are affected by the relativistic 
beaming effects. 

4.1. Source counts 

The MOJAVE sample is sampled from different high frequency 
radio samples and does not represent an independent survey. 
One needs to test its completeness which is the main crite- 
rion for deriving the true cosmological evolution and statisti- 
cal characteristics of relativistic jets of compact AGN. If the 
MOJAVE sample is a complete sample, then one should expect 
that the slopes of source counts of the MOJAVE sample and in- 
dependent high frequency surveys are aligned. At present, the 
only independent survey of ra dio sources at 15.2 G Hz is carried 
out using the Ryle Telescope fotvloret aljE oOR They identi- 
fied 66 sources between 0.2 Jy and 5Jy and fitted the differ- 
ential source count by a power law, n(S) oc S~ 2 . At 15.2 GHz 
and flux density range the main population of objects should 
consist of the beamed fl at-spectrum extragalactic radio sources 
dJackson & Wallll999l) . 

The normalized differential source count of AGN from the 
MOJAVE sample is presented in Fig. @] To take into account 
the different flux density limits of the sample, the northern 2n 
sky area is used to normalize the number count of AGN be- 
tween the flux density range from 1 .5 Jy to 2 Jy, and the num- 
ber count of AGN with S > 2 Jy is normalized to the complete 
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sky area of the sample, 2n(l + 0.342). The slope of the dif- 
ferential source counts (n(S) oc 5 , °Jy _I sr _1 ) is estimated for 
three sets of sources (Fig. 0): all AGN (a = -2.15 ± 0.11), 
quasars (a = -2.01 ± 0.09), and BL Lacs, galaxies and sources 
with no optical identification (a = -2.33 + 0.08). The MOJAVE 
and 15.2 GHz surveys overlap between the 1 .5 Jy and 5 Jy flux 
density range. The source counts of all MOJAVE sources (solid 
line) and quasars only (dashed line) have slopes which are char- 
acteristic of the complete survey at 15.2 GHz (a = -2) and of 
the bright end (> 1 Jy) of lower frequency surveys at 0. 15 GHz 

to 8.4 GHz (a 1.9) suggesting that the MOJAVE represents 

a flux-limited complete sample. The subsample of non-quasars 
(dot-dashed line) has a slightly steeper slope. It is most likely 
that this deviation arises because of a small number of bins (or 
sources) involved in the fitting. 

It is useful to compare the observed source count with the 
Euclidean source count which is given by N(> S) oc S~ 3 / 2 
for the unbeamed and uniformly distributed radio sources. It is 
not obvious whether this simple relation holds for the beamed 
sample of radio sources. Given the apparent radio luminosity 
function <J>(L, z), the number of radio sources with apparent flux 
> S per unit area of the sky is given by, 



l r L «« r 



,Z(S ' L) dV(z) 
dz 



(4) 



where dV is the comoving volume element. For the isotropic 
non-evolving Euclidean universe dV = 4jrR 2 dR and S = 
L/(4ttR 2 ), and the equation ® reduces to, 



N(> S)ozS 



-3/2 



r 



L V2 <$(L)dL. 



(5) 



The apparent luminosity function is the convolution of the in- 
trinsic luminosity function (iff(£)) an d probability density of a 
Doppler factor (P s (6)), 



L Up-i r 

®(L)dL = dh t(r(£)P s (L/£)d£, 

P J 



(6) 



w here the integration lim i ts are a com plex function of L and 
6 JUrrv & Padovanilll99lt lListerll2003l) . The integral over the 
apparent luminosity function is constant (eq. [5jl and hence the 
Euclidean differential source counts of a beamed population of 
radio sources, n(S) = dN(> S)/dS oc S~ 25 , is independent 
of the distributions of their intrinsic luminosities and Doppler 
factors. The slope of the differential number counts of the 
MOJAVE sources (a = -2.1) is flatter than the slope of the 
Euclidean source counts (a = -2.5) indicating that the radio lu- 
minosity function of compact AGN evolves with cosmic epoch. 

4.2. Radio luminosity function 

The observed radio luminosity of BL Lacs and quasars 
at 15 GHz is highly boosted and radiated predomi nantly 
by the a pproaching relativistic jet l lKeilermann et alJ 1 1998: 
IZensus et al.ll2002h . The relativistic Doppler beaming of the 
jet transforms the intrinsic radio luminosity function of com- 
pact quasars to the beamed (apparent) luminos ity function at 
the rest-frame of the observer. Earlier studies JUrrv & Shafer 
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Flux density S [Jy] 

Fig. 4. The differential source count of AGN sampled at 15 
GHz. Filled squared represent all radio sources, grey circles 
are quasars and filled diamonds are BL Lacs, galaxies and un- 
known type of sources. The length of horizontal bars marks the 
flux density bin size. 



Il984l lUrrv & Padovanilll99lt iListerl 120031) provided analyti- 
cal predictions for the apparent luminosity function of jets ran- 
domly oriented on the sky. The single and double power-law in- 
trinsic luminosity functions were modelled assuming a single- 
valued Lorentz factor for the jet. IListerl y003) extended the 
later work for the variety of power-law distributions of Lorentz 
factors, P y oc y (y] < y < -y 2 ), which was found to pro- 
yide a good fit to t he distribution of apparent speeds of AGN 
dLister & Marscherl 19971) . The modelled apparent and intrinsic 
luminosity functions have several important break points and 
relations. We will use these to constrain the intrinsic luminos- 
ity function of the jets. The evolution of the apparent luminos- 
ity function involves both evolution of the intrinsic luminosity 
function and evolution of the Lorentz factor (or the change of 
the Doppler factor) of the jet with redshift. The latter can be 
examined by considering the apparent speed-redshift relation 
plane which also allows the important characteristics of the rel- 
ativistic beaming of flat-spectrum AGN to be determined. 

4.2.1 . The apparent speed - redshift plane 

The apparent speed j6 a (in units of the speed of the light, c) of 
jets measured for 87 AGN is plotted against their redshifts in 
Fig-0 The mean apparent speed varies significantly depending 
on type of AGN: (1.9 ± 1 .0)c for radio galaxies, (5.9 + l.l)cfor 
BL Lacs and (9.4 ± 1.0)c for quasars. Most of quasars and BL 
Lacs have superluminal jets (j8 a > 1) indicating that their radio 
emission is highly beamed. This puts a limit on the minimum 
Doppler factor, <5 m in = 5(y\ , 9 C ) — 1 (L <: £.), where 9 C is the 
maximum (critical) viewing angle of the jet having a m inimum 
Loren tz factor j\. Assuming that 9 C w 45° for quasars ( B arthell 
1989), we calculate the lower limit on the Lorentz factors, j\ ~ 
3 ( from eq.B). 

IVermeulen & Cohenl d 19941) showed that for the 
orientation-biased sample with a single Lorentz factor 
distribution the value of the apparent speed is close to the 
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Redshift z 



Fig. 5. Apparent transverse speed of the jet as a function of red- 
shift for AGN with known values of y8 a : 66 quasars (filled cir- 
cles), 15 BL Lacs (gray diamonds) and 6 radio galaxies (filled 
triangles). The gray squares are the mean apparent speed of 
quasars estimated in six redshift bins (binning ranges are the 
same as in Fig. [6}. The dotted line shows the division of quasars 
in two equal subsamples with low (J3 d < 6) and high apparent 
speeds Q3 d > 6). 1 <x error are presented for all data. 
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Fig. 6. The cumulative redshift distribution of compact quasars. 
The z-distribution of 66 quasars with known values of apparent 
speed (grey area) is superimposed with the distribution of 28 
quasars with unknown /3 a (clear area). 

maximum apparent speed of the jet y8 a w y8 ajtnax = yjy 1 — 1 as y 
which is true for BL Lac objects and quasars with y 2 » 1, 
and (5 max ~ 2y2 = 2j8 amax . The maximum apparent speed 
and Lorentz factor of the jet is 72 ~ Ai.max ~ 30 (Fig. 13 
and the £ max » 60. The maxumum Doppler factor will be 
reached at the viewing angle * degrees. It is extremely 
unlikely to detect a source with £ max = 60 in a small sample 
such as the MOJAVE. 5 ~ j6 a , max at the optimum angle of 
the jet, sin 6 = l/y, and the value of 5 ~ 30 is a realistic 
upper limit to the Doppler factor of the jets of our sample. 
To calculate the range of the intrinsic luminosities of jets we 
use the minimum and maximum radio luminosities of quasars 
L min = 2.3 x 10 26 WHz 1 (atz ~ 0.2) andL max = 10 29 WHz 1 
(at z ~ 2). Assuming that 8 ~ /? a , max ~ 10 at low and high red- 



shifts (see Fig.© we recover L\ = eT 2 L min ~ 2.3xl0 24 WHz 1 
and£ 2 = <T 2 L max * lO^WHz- 1 . 

There is a lack of high /? a at low redshifts, z < 1 (Fig.|5|i. For 
the j6 a - L relation plane, Cohen et al. (in preparation) discuss 
the selection effects which might result in the lack of sources 
with large /3 a at low luminosities. They argue that both effects, 
the minimum sampling interval of the MOJAVE sample and 
the determination of speed of the jet when the frequency of 
ejection of moving components is greater than the frequency 
of observations, can not account for the lack of sources thus 
leading to the conclusion that the increase of apparent speeds 
with increasing luminosity is real. This is true also for the j3 d -z 
relation plane because of positive correlation between L and z 
in the flux-limited MOJAVE sample. Another potential selec- 
tion effect at low redshifts is that we miss sources with un- 
known p. d which may populate the left-top region in Fig. [5] The 
quasars with unknown apparent speeds (Fig. |6jl will populate 
the redshifts z > 0.5 thus the increase of y6 a .max to high red- 
shifts (z ~ 1) will not be affected. This trend is also evident for 
the mean apparent speeds (gray squares in Fig. [5} estimated in 
different redshift bins. There is a gradual decrease of the mean 
apparent speed from z ~ 1 down to z ~ 2.5. Were high speed 
jets (y3 a > 30) to exist at high redshifts one should observe 
them because their radio emission would be highly boosted. It 
is unlikely that the high-speed quasars are among those with 
unknown /3 d , because no systematic selection of measured ap- 
parent speeds was introduced. For the same reason, the aver- 
aged behaviour of the apparent speed can not be affected by 
unknown apparent speeds. Moreover, the highest redshift bin, 
which includes six quasars with measured j0 a , has the lowest 
value of (J3 a ) = 5.2 ± 1 .7 which supports the idea that the appar- 
ent speeds (and hence the Lorentz factors) of the jets of quasars 
evolve with cosmic epoch: the mean apparent speed increases 
from z ~ 2.5 to z ~ 1 by a factor 2.5 and then falls by factor of 
two at the present epoch. 

Note that a peak in the (fi d ) - z relation plane is mainly due 
to nine high speed sources with /? a > 20, and that inclusion 
or exclusion of these quasars has no significant effect on the 
evolution of the luminosity function of quasars examined in 
the next sections. 

4.2.2. Evolution of the luminosity function of quasars 

We chose the population of quasars for a study of the lumi- 
nosity function because they comprise a large homogeneous 
sample of beamed and optically identified sources having ap- 
pare nt speeds m easured for ~70% of the jets. We used the 1 / V a 
(Schmidt 1969) method to construct the apparent luminosity 
function. The MOJAVE sample occupies different sky areas de- 
pending upon the flux density limits: 1.5 Jy for northern hemi- 
sphere and 2 Jy for southern hemisphere. In thi s case, the V a is 
the vo lume available in the combined sample dAvni & Bahcalll 
1980). Note that all statistical tests and methods throughout the 
paper were applied to the combined sample. The apparent lu- 
minosity function is presented in Fig.0for three redshift bins, 
< z < 0.7 (27 quasars), 0.7 < z < 1.4 (35) and z > 1.4 (32). 
The form of the luminosity function is similar at all redshifts: 
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27 27.5 

Log L 15 / W Hz" 1 

Fig. 7. The radio luminosity function of 94 quasars is presented 
for three redshift bins z < 0.7 (left curve), 0.7 < z < 1.4 
(middle curve) and z > 1.4 (right curve). The Poisson error 
(1 cr) acossiated with the luminosity function is given (cr 2 = 
(ALogiL^Y^HVl). 



it is fiat at lower luminosities (A ~ 0, for <$>(L) = L~ A ) and be- 
comes steeper A ~ (1.9 to 2.5) at higher luminosities. There is 
marginal evidence that the high luminosity end of the luminos- 
ity function steepens with increasing redshift. We now discuss 
the form of the intrinsic luminosity function and distribution 
of Lorentz factors which may reproduce the observed luminos- 
ity function. We consider a single (steep) and double (fiat and 
steep) power-law intrinsic luminosity functions (0(X) = £T B , 
where £\ < £ < £?.) a nd distribution of Lorentz factors with 
yi < y < 22 dUrrv & Padovanill99ltlListerl2003l) . 

lUrrv & Padovanil Jl99ll) showed that if the slope of the ap- 
parent luminosity function is steeper than -(p + l)/p = -1.5 
(p = 2) then the intrinsic and apparent luminosity functions 



have the same slope at high luminosities L > L\ — & p mwi £.\ 
(here we use the notation of Urry & Padovani 1991). At high 
luminosities the slope of the apparent luminosity function (and 
henc e the slope of the i ntrinsic luminosity function) is ~ 2.5. 
IUrrv&Padorarllll99lh claimed that such a single and steep 
power-law intrinsic luminosity function will be transformed to 
the apparent luminosity function which, at low luminosities 
(L < L4), will have a slope — (p + l)/p = -1.5 independent 
of the Lorentz factor distribution. A single power-law parent 
luminosity function can not produce the fiat part (B ~ 0) of 
the apparent luminosity function observed at low luminosities. 
iListerl J20 03) showed that a fiat or inverted slopes may be pro- 
duced with the steep power-law slope (B - 2) of the intrinsic 
luminosity function, 1 < y < 30 and steep (a < 0) distribu- 
tion of the Lorentz factors. This is true for orientation-unbiased 
samples populated by slow jets and apparently it is not appli- 
cable to the flux-limited MOJAVE sample which is dominated 
by relativistic jets and high 5 > 1 . 

A complex (fiat and steep) power-law intrinsic luminosity 
function and single power-law distribution of Lorentz factors 
of the jets are needed to describe the observed flatness of the 
luminosity function of beamed quasars. At low luminosities, 
the fiat slopes of the intrinsic and apparent luminosity func- 



tions should coincide llUrrv&Padovanilll99lh while the steep 
slopes should match at high luminosities. The overall intrin- 
sic luminosity function can be crudely represented by a double 
power-law with </>(£.) ~ const and (p(X.) ~ XT 25 at low and 
high luminosities respectively. The broken power-law form is 
generally consistent with the luminosity f unctions of quasars 
selected from unbeamed surveys: optical (Boyle et al. 1988) 
and low radio frequency samples ( Willot t et alJll998L EoOl). 
This supports the idea that fiat-spectrum quasars are the parent 
population of FRII radio sources beamed towards the observer 
jjackson & Walll999h . 

The shift of the apparent luminosity function (Fig.0 may 
happen as a result of the evolution of intrinsic luminosity func- 
tion of quasars and/or evolution of the Lorentz factor (Doppler 
factor) of the jet with cosmic time. We now consider three 
cases: a pure evolution of the Lorentz factors of the jets, a 
pure evolution of the intrinsic luminosity function, and evo- 
lution of both. The mean Lorentz factors in three redshift bins 
(Fig. 01 are (8.3 + 1.4) at low-z, (11.6 ± 2.2) at intermediate- 
z and (8.4 + 1.8) at high-z. In the case of an evolution of the 
jet speeds, the unchanged intrinsic luminosity function will be 
transformed to the apparent luminosity function. The observed 
shift of the luminosity function at high redshifts requires the 
increase of the beaming factor (or the Lorentz factor) with in- 
creasing redshift. The mean Lorentz factor increases from low- 
z (8.3 ± 1.4) to intermediate-z (11.6 ± 2.2) and goes down at 
high-z (8.4 ± 1.8). This is in contradiction with the expected 
increase of the beaming factor at high redshifts. 

In the case of a pure luminosity evolution of the intrin- 
sic luminosity function one should expect that the luminosity 
curves shift to high luminosities (Fig. [7) at high redshifts (sim- 
ilar to the shift of the luminosity func tion of opti cally selected 
quasars, e.g. Boyle et al. 1988~). IWall et all d2005l) noted that a 
contamination of the low luminosity bin by unbeamed sources 
can mimic a luminosity evolution. Our sample includes only 
the beamed quasars and therefore the observed shift is likely 
to be produced by simple luminosity evolution. An important 
question is whether the evolution of the beaming effects have 
a contribution to the shift of the luminosity function. Although 
there are indications that the Lorentz factor evolves with cos- 
mic time (Fig. [5} it has no significant effect on the evolution 
of the apparent luminosity function. To understand and recover 
fully the effects of the relativistic beaming on the intrinsic lumi- 
nosity function one needs to compare the observed distributions 
of jet parameters (and a complete data set of apparent speeds) 
with the distributions generated by means of Monte-Carlo sim- 
ulations (Lister & Marscher 1997; Lister et al., in preparation). 

4.3. V/V max evolution 



Here , we use the model-independent V/ V m ax statistics dSch midt 
1969) to test a null hypothesis of a uniform distribution of com- 
pact AGN in space and to study their cosmological evolution. 
This test is applicable only for complete flux-limited samples, 
where V is the volume of space enclosed by the redshift of 
the AGN having a certain flux density, and V max is the maxi- 
mum volume of space within which this source could be ob- 
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Table 1. V e /V a test for 1 16 AGN. 



ID 


Num 


(VJV,) 


Standard 


K-S test 






N 




Error a 


(%) 


Prob. 


>° 
i 


All 


116 


0.581 


0.025 


98.83 


0.012 


>° 

\ *> 


Quasars 


94 


0.589 


0.027 


99.02 


0.009 


n. ° 
> 


BL Lacs 


15 


0.617 


0.070 


52 


0.48 


i 


Galaxies 


7 


0.375 


0.102 


68 


0.32 





1 The standard error is estimated from SE = ^j(x 2 ) - (x) 2 / yfFJ, 
where x = V e /V a and N is the number of sources. 




served and still be included in the complete sample. If sources 
have a uniform distribution in space, then the averaged ratio 
of V/ V m ax is constant and equal to 0.5. The difference between 
the expected and observed values of V/V max allows the signifi- 
cance of a deviation of the real distribution from uniform to be 
estimated. This test has been applied for different radio sam- 
ples to show the presence of a strong posi tive evolution up to 
z ~ 2 with (VIV m ^) ~ 0.6 to 0.7 (e.g JLongair & ScheueJ 
1 970; IWillsfcLvndsll978t IPeacock et all 198 



W e employed the generalized V/ V max test I Avni & Bahcalll 
1980) which is derived for combined samples with different 
flux-density limits and sky areas. The new variables are V e and 
V a , the volume enclosed by th e source and the volum e avail- 
able in the combined sample. IXvni & Bahcalll Jl980l) showed 
that the formal statistical properties of V e /V a are the same as 

those Of V/Vmax- 

For 94 quasars, <V e /V a > = 0.589 ± 0.027, the distribution 
is not uniform at the 99.3 % confidence level (see Table 2), 
and it is biased towards large values which indicates that com- 
pact quasars, and hence the powerful jets, were more populous 
at high redshifts. The BL Lacs show a similar trend, but this 
is not statistically significant because of the small number of 
sources. The distribution seems to be uniform for galaxies, i.e. 
with no evolution so far. The plausible explanation is that all 
eight radio galaxies occupy the low-redshift region where the 
density/luminosity evolution is negligible, but better statistics 
would be needed to confirm this result. 

To investigate the dependence of radio luminosity on the 
Ve/Va statistics, we divided the sample of quasars in two 
equal subsamples by absolute luminosity at 15 GHz, L\$ = 
10 27 " 9 WHz -1 . Strong sources have redshifts greater than 0.5 
and weak sources have z < 1.5 (Fig.Q. For 46 low-luminosity 
quasars we find that <V e /V a } = 0.658 ± 0.036 with confi- 
dence level of 99.96 %, indicative that the distribution of V e /V a 
is biased towards large values, whilst for 48 strong sources, 
<Ve/Va) = 0.53 ± 0.04 (41 %), no significant deviation from a 
uniform distribution was found. The K-S test rejected at 96 % 
confidence level the null hypothesis that these distributions are 
drawn from the same parent population. The Student t test re- 
jected the hypothesis of equal mean values of V e /V a for low- 
and high-luminosity quasars at high significance level, 0.011 
(98.88%). This suggests that the cosmological evolution of 
quasars depends on radio luminosity. The evolution of high- 



Fig. 8. The banded V/V max statistic <(V e - V )/(V a - V )> against 
Zq, the lower limit in the redshift bins. This plot is presented 
numerically in the Table 3. lcr standard error (see Table 3) is 
presented. 



redshift quasars can be masked by the strong positive evolu- 
tion at low redshifts. To investigate the evolutionary behaviour 
of quasars at high redsh ifts we used the banded version of the 
V/V max test JOsmer & Smith! 197(11 Avni A Shilleill983h which 
allows one to mask out evolutionary effects at high redshifts. It 
calculates the mean of a new statistic, (V e - Vo)/(V a - Vo), for 
the samples restricted by different values of zo (z > Zo), where 
Vo is the volume enclosed by a redshift zo- 

The results of a banded V/Vmax test are shown in Figs. 181101 
For 94 quasars, a weak positive evolution is seen at low red- 
shifts, z < 1, with values greater than the 0.5 no-evolution line 
(Fig. |8j. The positive evolution is significant for quasars with 
Z < 0.2 (see Table 3). The values of (V e /V a ) decrease gradually 
to 0.5, and, within the errors, they remain at no-evolution level 
out to z ~ 1.6. A negative evolution is becoming significant for 
quasars with z > 1.7 (Table 3), and it remains significant out 
to z ~ 2.5 at which only few sources can be found. A signif- 
icant redshift cut-off indicates that the population of beamed 
quasars having powerful jets undergo a diminution already at 
Z ~ 1 .7. Slightly high redshif t cut-off (z > 1 .8) is fou nd for faint 
quasars sampled at 1 .4 GHz t Cirasu oio et aljEo05l) and for the 
combined flat-spectrum quasars at 2.7 GHz (z > 2, Dunlop & 
Peacock 1990). 

Figure [9] shows the results of a banded test for low- and 
high-luminosity quasars separated by luminosity at L15 ghz = 
10 27 ' 9 WHz -1 . There is significant positive evolution for low- 
luminosity quasars out to z ~ 0.5, and, again, a strong evi- 
dence for a high-redshift decline in the comoving space density 
of high-luminosity quasars with z > 1.7. The luminosity ap- 
pears to play an important role in the evolution of flat-spectrum 
quasars. 

We find that the evolution of flat-spectrum quasars depends 
strongly on the speed of the jet. Relatively slow jets (J3 a < 6, 
see Fig. [5} do not show significant evolution up to z ~ 2.5 
(Fig. llOt . whilst the population of fast jets grows rapidly be- 
tween z ~ 2.5 and z ~ 1.7 and remains at no evolution level at 
small redshifts. As discussed in Sect. 4.2.1, the missed quasars 



8 



T. G. Arshakian et al.: Cosmological evolution of compact AGN 



Table 2. Tabulated (Ve/V^} results for 94 flat-spectrum quasars. 



Redshift 
range 



Number 

N 



(VJV. d )± SE a 



K-S test 
Prob. 



>0.1 


94 


0.592 


± 0.027 


0.008 


>0.3 


90 


0.583 


± 0.029 


0.047 


>0.5 


84 


0.557 


±0.031 


0.280 


>0.7 


67 


0.542 


± 0.036 


0.361 


>0.9 


56 


0.547 


± 0.041 


0.295 


>1.1 


43 


0.543 


± 0.043 


0.383 


>1.3 


33 


0.494 


± 0.046 


0.937 


>1.5 


25 


0.462 


± 0.048 


0.352 


>1.7 


20 


0.380 


± 0.053 


0.021 


>1.9 


13 


0.319 


±0.071 


0.009 


>2.1 


8 


0.301 


±0.105 


0.012 


>2.3 


5 


0.303 


±0.177 


0.114 



1 The standard error: SE = V<* 2 > - <*> 2 / VlV, where 
x = (V e - Vo)/(V a - Vo) and N is the number of sources withz > zo- 



Low L,.< 10" WH! 



□ Highl,.> 10" WHz 



0.5 1 1.5 2 2.5 0.5 

Z 



1 1.5 2 2.5 

Z 



Fig. 9. The banded V/V max statistic versus zo for 46 low- 
luminosity (left panel) and 48 high-luminosity (right panel) 
quasars. 



with no y6 a values should populate randomly the samples of 
high- and low-speed quasars and hence it is unlikely that their 
exclusion will introduce the observed differential evolution of 
quasars with high and low speeds. For a given apparent lumi- 
nosity the population of slow jets should have brighter intrinsic 
luminosities than those of slow jets. This hints that there should 
be two distinct populations of quasars, one with fast jets and 
low intrinsic luminosities which forms at later epochs and un- 
dergoes a strong evolution up to z ~ 1 and one with slow jets 
and high intrinsic luminosities formed at much earlier cosmic 
epochs and evolving slowly to present time. A dual population 
model consisting of an old population of slowly evolving low 
speed jets and young population of strongly evolving fast jets 
at z > 1 may produce an apparent 'negative evolution' of yS a 
seen at z > 1 (Fig. 0. In this scenario, the positive evolution 
of the apparent speed at z ^ 1 (Fig. [3 is mainly due to young 
population of high speed quasars. 



dio luminosity function is fiat (B ~ 0) at low luminosities 
and becomes steeper (B ~ -2.5) at high luminosities; (f) the 
weak positive evolution of quasars is evident at low redshifts, 
z < 0.5, and a statistically significant redshift cut-off is present 
at z > 1.7; (g) the cosmological evolution is different (1) 
for high- (L15 > 10 27 9 WHz _1 ) and low-luminosity quasars 
(Z45 < 10 27 9 WHz _1 ) suggesting that their evolution depends 
strongly on luminosity, and (2) for fast and slow jets, which in- 
dicates the existence of two distinct populations of quasars one 
with slow jets (/? a < 6) and high intrinsic luminosities which 
does not show a significant evolution with redshift and one with 
j8 a > 6 and low intrinsic luminosities having a pronounced red- 
shift cut-off at z > 1.7. 

The differential evolution of the population of fast and slow 
jets at z > 1.7 may explain inconsistent results obtained for 



5. Conclusions and discussions 



The principal results of statistical analysis of bright fiat- 
spectrum compact AGN from the flux-limited MOJAVE sam- 
ple are: (a) the distribution of 133 AGN is uniform on the sky; 
(b) the slope of the source counts of radio sources is aligned 
with the slopes of source counts of independent high radio 
frequency surveys suggesting that the MOJAVE is a complete 
flux-limited sample; (c) Lorentz factors of the jets of quasars 
lie in the range from ~ 3 to ~ 30 and Doppler factors of the 
jets are distributed between £ m ; n > 1 and £ max ^ 60 with the 
realistic upper limit on 6 max « 30 for the MOJAVE sample; 
(d) the mean apparent transverse speed of jets of quasars (or 
the mean Lorentz factor) increases by a factor of two (from 
6.5 to 14) with increasing redshift up to z ~ 1 and then grad- 
ually decreases at high redshifts by a factor of 2.5; (e) the in- 
trinsic luminosities of compact jets are distributed in the range 
from ~ 10 24 WHz _1 to ~ 10 27 WHz 1 and their intrinsic ra- 



A 
>° 

>" 



> 

V 




Fig. 10. The banded V/V max statistic versus zo for 32 low-j6 a 
(left panel) and 34 high-/? a (right panel) quasars with known 
apparent speeds. 
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a redshift c ut-off of quasars sampled at high and low radio 
frequencies JPunlop & Peacock! 1990lrWillott et all200ll) with 
former one providing no evidence for any decline in the comov- 
ing space density. If slow jets reside in high luminosity FRII 
radio sources and fast jets produce faint FRII structures, which 
are too weak to be detected at 15 1 MHz, then the low-frequency 
sample will include the population of quasars with low jet 
speeds which do not show significant evolution up to z ~ 2.5 
(Fig. I10> . The population of fast jets, however, will be in- 
cluded in the high frequency samples because of highly beamed 
emission of the jet. This naturally explains the reds hift cut-off 
of quasars sampled at 1 .4 GHz, 2.1 GHz JPunlop & Peacock! 
ll990HCirasuoloetalJl2005UWalletalJl2005t) and significant 
redshift cut-off at z > 1.7 found for our sample ("Figs. l9TTUI >. 
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